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ABSTRACT 

We have investigated the critical conditions required for a steady propeller effect for 
magnetized neutron stars with optically thick, geometrically thin accretion disks. We 
have shown through simple analytical calculations that a steady-state propeller mech¬ 
anism cannot be sustained at an inner disk radius where the viscous and magnetic 
stresses are balanced. The radius calculated by equating these stresses is usually found 
to be close to the conventional Alfven radius for spherical accretion, ta- Our results 
show that: (1) a steady propeller phase can be established with a maximum inner disk 
radius that is at least ^15 times smaller than ta depending on the mass-flow rate of 
the disk, rotational period and strength of the magnetic dipole field of the star, (2) 
the critical accretion rate corresponding to the accretion-propeller transition is orders 
of magnitude lower than the rate estimated by equating ta to the co-rotation radius. 
Our results are consistent with the properties of the transitional millisecond pulsars 
which show transitions between the accretion powered X-ray pulsar and the rotational 
powered radio pulsar states. 

Key words: pulsars: individual (PSR J1023-I-0038, XSS J12270—4859) - accretion 
- accretion disks 


1 INTRODUCTION 

Rotational properties of a magnetized neutron star evolving 
with an accretion disk are governed by the interaction be¬ 
tween the disk and the magnetic dipole field. The star could 
spin-up or spin-down depending on the mass-fiow rate, Min, 
of the disk, strength of the magnetic dipole field, B, on the 
surface of the star and the details of the disk-magnetosphere 
interaction. The inner disk radius, rin, is usually estimated 
by equating the viscous and the magnetic stresses and found 
to be close to the conventional Alfven radius, ta, which is 
calculated for spherical accretion by equating the ram pres¬ 
sure of matter to the magnetic pressure (Lamb et al. 1973, 
Davidson & Ostriker 1973). When rin is less than the co¬ 
rotation radius, rco, at which the field lines co-rotate with 
the disk matter, the system is in the spin-up phase, and the 
mass from the inner-disk can flow onto the star along the 
field lines. When rin is greater than rco, the star is in the 
spin-down state, which is also known as the ‘propeller’ phase 
(Illarionov & Sunyaev 1975), since in this regime, it is usually 
assumed that the inflowing disk matter is expelled from the 
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system by the dipole field lines rotating with speeds greater 
than the local escape speed. 

Observations of three recently discovered transitional 
millisecond pulsars (MSPs) which show transitions between 
the accretion powered X-ray pulsar and the rotational pow¬ 
ered radio pulsar states provide an excellent laboratory to 
study the physics of the disk-magnetosphere interaction (see 
e.g. Linares 2014 for a recent review). X-ray pulsations from 
the transitional MSPs XSS J12270—4859 (Papitto et al. 
2015) and PSR J1023+0038 (Archibald et al. 2015) were 
detected during sub-luminous accretion states with X-ray 
luminosities of a few 10^® — 10®^ erg s“^ which correspond 
to the accretion rates of ~ 2 — 5 x 10^® g s“^. The most 
plausible origin of these X-ray pulses seems to be accretion 
onto the poles of the star (Papitto et al. 2015, Archibald et 
al. 2015). For B — 2 y. 10® G, a typical field estimated for 
the millisecond pulsars with periods of a few ms (Lyne & 
Graham-Smith 2007; see also Section 2.3), ta is inferred to 
be greater than even the light cylinder radius, rue = c/D* 
where c is the speed of light, and fl* is the angular fre¬ 
quency of the neutron star. This is in sharp contrast to the 
expected condition, rin — rco, for accretion onto the star. A 
large fraction of the inflowing disk matter might be thrown 
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out of the system, while the observed X-ray luminosity is 
powered by the remaining few per cent of the disk mass- 
flow accreted onto the neutron star. In this case, since rin 
could come close to rco accretion is possible, but the X-ray 
luminosity of the disk would be much higher than the ob¬ 
served luminosity. Is it possible to overcome this difficulty 
by invoking radiatively inefficient flow which could prevail 
when the inner disk becomes optically thin (Papitto et al 
2015). This does not seem to be likely, because the disk be¬ 
comes optically thin only when Min decreases below ~ 10^^ 
g s“^ (Done et al. 2007). The inner disk is expected to be 
optically thick at Min — 5 x 10^® g as inferred from 
the observed X-ray luminosity assuming no outflow. With a 
much higher Min needed for rin — rco, the inner disk would 
have an even higher optical thickness. So, one cannot invoke 
an optically thin disk to avoid the problem of unobserved 
disk luminosities if rin — rco. The inferred rin — rA in tran¬ 
sitional MSPs is in fact found to be even greater than rue- 
So, how can these transitional MSPs accrete matter from 
the disk? 

In the spin-down phase, if the inflowing disk matter 
cannot be repelled from the system efficiently, the resul¬ 
tant pile up at the inner disk increases the density and 
the temperature in the disk, rin propagates inwards, and 
eventually reaches the co-rotation radius. Resultant onset 
of accretion could evacuate the inner disk, then rin moves 
outwards stopping the accretion. Subsequent refilling of the 
inner disk resuming the accretion could lead to a cyclic evo¬ 
lution of the inner disk even when the mass transfer from 
the outer disk is constant. This idea was first proposed by 
Sunyaev & Shakura (1977) and Spruit & Taam (1993), and 
recently developed through a series of papers by D’Angelo 
& Spruit (2010; 2011; 2012). In these works, it is assumed 
that the propeller mechanism is not efficient even for very 
low disk mass-flow rates. With this assumption, numerical 
simulations show that, for a large range of initial conditions 
including those with very low disk mass-flow rates, sources 
evolve in the long-term into the so-called trapped states, 
with rin — rco. These states could be persistent or the inner 
disk could go through cyclic accretion episodes on timescales 
between the dynamical and viscous timescales depending on 
the average mass-flow rate from the outer to the inner disk 
(D’Angelo & Spruit 2012). 

The accretion and rotation powered states of the tran¬ 
sitional MSPs do not seem to be the results of the cyclic 
processes described by D’Angelo & Spruit (2012). The prop¬ 
erties of these sources imply that there is another mechanism 
that impedes accretion below a critical disk mass-flow rate. 
The radio pulses are always observed below a certain X-ray 
luminosity less than a few 10®^ erg s“^, while Lx > 10^® 
erg s“^ when the X-ray pulsations are observed in the sub- 
luminous X-ray states. The sources remain in one or the 
other state on timescales (months) much longer than the dy¬ 
namical or viscous timescales of the inner disk. This might 
indicate that the propeller mechanism is at work below a 
certain accretion rate. Flat radio spectra observed in the 
sub-luminous X-ray states could be a sign of the mass out¬ 
flows (Papitto et al. 2015, Linares et al. 2014). The radio 
and X-ray pulses have not been observed in the same epoch, 
or at different epochs with the same X-ray luminosity. This 
strongly indicates that the radio pulses are switched on when 
the accretion is impeded, while the X-ray pulses are pro¬ 


duced as a result of mass-flow onto the star, which switches 
off the pulsed radio emission. 

In the disk-magnetosphere interaction model employed 
by D’Angelo & Spruit (2010; 2011; 2012), the magnetosphere 
is formed by the closed field lines that extend up to a ra¬ 
dius, rm. The interaction between the disk and the dipole 
field takes place inside this boundary region between Vm and 
Tm -I- Ar, while the disk is decoupled from the field beyond 
the interaction region. Inside the boundary, the closed field 
lines cannot slip through the disk, since the magnetic diffu¬ 
sion timescale, which is comparable to the viscous timescale 
(Fromang & Stone 2009), is orders of magnitude longer than 
the disk-field interaction timescale. The field lines interact¬ 
ing with the inner disk inflate and open up. This disk-field 
interaction model with a narrow boundary was proposed by 
Lovelace, Romanova & Bisnovatyi-Kogan (1995), and sup¬ 
ported by the numerical simulations (Hayashi, Shibata, & 
Matsumoto 1996, Goodson, Winglee, & Boehm 1997, Miller 
& Stone 1997). We will examine the propeller effect in the 
context of this disk-magnetosphere interaction model which 
was developed further to study the disk-field interaction and 
the outflows in the propeller phase (Lovelace et al. 1999, 
Ustyugova et al. 2006). We will use the term ‘propeller 
state’ to indicate the phase in which the accretion onto the 
star is not allowed and the disk matter arriving at the disk- 
magnetosphere boundary is thrown out of the system with 
speeds greater than the escape speed. 

In the present work, we investigate the physical condi¬ 
tions required for a steady or long-lasting propeller phase. 
In Section 2, we derive the maximum inner disk radius in 
a propeller regime and the critical accretion rate for the 
accretion-propeller transition. We discuss our results with 
comparisons to the characteristic properties of the transi¬ 
tional MSPs, and summarise our conclusions in Section 3. 


2 PROPELLER PHASE 

2.1 Critical condition for the propeller mechanism 

In the case of spherical mass-flow onto a neutron star, the 
matter falls with the free-fall speed, vg. The spherical mass- 
flow rate can be written as 

Mg = 47rr^ps(r) Ufi(r) (1) 

where ps{r) is the density of the matter at radial distance r. 
Magnetic dipole field of the star can stop the mass-flow at 
the radius where the pressure of the falling matter becomes 
equal to the magnetic pressure 

^ = ps{r) vg (2) 

where B = p/E, and p is the magnetic dipole moment of 
the star, equation (2) is satisfied at the Alfven radius 

rA ~ (GM)-^/'^ (3) 

(Lamb et al. 1973, Davidson & Ostriker 1973) where G is 
the gravitational constant, and M is the mass of the neutron 
star. The inner disk radius for a geometrically thin accre¬ 
tion disk is is usually assumed to be the radius at which the 
viscous and the magnetic stresses are balanced, and calcula¬ 
tions give inner disk radii very close to nominal Alfven radius 
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va given by equation (3) within a factor of ~ 2 (Spruit and 
Taam 1993). 

In the present work, we adopt the basic principles of the 
disk-magnetosphere interaction model proposed by Lovelace 
et al. (1995), and developed later by Lovelace et al. (1999) 
and Ustyugova et al. (2006) to study the propeller regime of 
magnetized stars in more detail. In this model, the radius of 
the axisymmetric magnetosphere, rm, is defined as the ra¬ 
dius of the region of closed field lines. Inside Cm, all matter 
and the field lines rotate together with the same angular ve¬ 
locity. The disk interacts with the large-scale field mainly in 
a narrow boundary region from r^n to rmT Ar. In the bound¬ 
ary region, azimuthal component of the magnetic field, B^, 
is generated through interaction with the inner disk, and 
is responsible for the torque acting on the disk. Since the 
diffusion timescale of the field lines are much longer than 
the interaction timescale of the disk matter and the field, 
grows to a maximum, with B^ « B^, and subsequently 
the field lines inflate and open up (Aly 1985, Lovelace et al. 
1995, Hayashi et al. 1996, Miller & Stone 1997, Uzdensky 
et al. 2002, Uzdensky 2004). The matter accelerated by the 
magnetic torque leaves the system along the open field lines. 
Reconnection of the field lines passing through the boundary 
layer leads to a new cycle of interaction, inflation, opening 
and reconnection. The outflow of the matter along the field 
lines is a complicated process and not well understood (see 
e.g. Uzdensky 2004 for a review). Here, we will simply as¬ 
sume that the acceleration of matter, within the interaction 
timescale, to the speed of the field lines rotating faster than 
the escape speed inside the boundary region is the minimum 
critical condition for an efficient propeller effect. 

The timescale of interaction between the field and the 
disk matter tint — — flK(nn)|~^ where fl* is the angular 

frequency of the neutron star, and IlK(rin) is the Keplerian 
angular speed at r = rin. If the mass inflow rate into the 
boundary, Min, is constant and the reconnections occur suf¬ 
ficiently rapidly leading to an average mass-outflow with a 
rate Mout = Min, the system could be well approximated 
by a steady-state model, neglecting the small-scale rapid 
variations in the boundary. Indeed, the simulations indicate 
that the inflation-reconnection events occur on the dynami¬ 
cal timescale at the inner disk, ~ 27r/flK(rm) (Lovelace et al. 
1999, Ustyugova et al. 2006). In the propeller phase, when 
the outflowing matter does not return back to the disk, and 
the outflows from outside the boundary are negligible, the 
disk outside tin = rm + Ar is expected to have the prop¬ 
erties of a standard thin disk. For a narrow Ar that is a 
small fraction of rm, the inner radius of the unperturbed 
disk rin = rm. 

There are a few important points that should be taken 
into account in the calculations of the the inner disk ra¬ 
dius in the propeller phase to remain self-consistent with 
the assumed nature and timescales of the mass outflows. 
The matter inside the boundary is accelerated from the Ke¬ 
pler speed to the speed of the field lines moving with speeds 
greater than the escape speed on the interaction timescale, 
that is, the speeds do not increase gradually from vk at rin 
to the speed of the field lines at rm along many orbits. In¬ 
stead the matter should be accelerated rapidly, because the 
interaction timescale corresponds to a small fraction of the 
local Keplerian period. Depending on the efficiency of recon¬ 
nections of the field lines, a given piece of matter couples to 


the reconnected field lines at a point inside the boundary, ac¬ 
celerated to high speeds in the azimuthal direction, and ex¬ 
pelled form the system along the open field lines. As we show 
in Section 2.2, this radius cannot be obtained from the bal¬ 
ance between the magnetic and the viscous stresses. Instead, 
the magnetic stresses inside the boundary must satisfy the 
condition imposed by the interaction timescale. There is a 
certain amount of angular momentum that should be trans¬ 
ferred to the matter to increase its speed above the critical 
level before the currently interacting closed field lines open 
up. The minimum magnetic torque required for this angu¬ 
lar momentum transfer depends on how fast this should be 
done, that is, on the interaction timescale. 

2.2 Escape radius 

To estimate the maximum inner disk radius in the propeller 
phase, we take B^ ~ B^, and B^B^ « B^ = [j? jA which 
gives the maximum attainable magnetic torque inside the 
boundary. With this assumption, magnetic torque per unit 
radial distance of the disk becomes B'^A' (Ghosh & Lamb 
1979). Cross sectional area of the disk perpendicular to the 
radial direction is 27rr 2/i where H = 2h is the full thickness 
of the disk. The magnetic torque per unit volume of the 
disk can be written as re = i?^r/(47r/i) (e.g. Rappaport 
et al. 2004). We neglect any heating effect by the propeller 
mechanism inside the boundary. In the steady state, mass- 
flow rate of the geometrically thin disk at a radial distance 
r can be written as 

M = 2nr E(r) Vr = 27rr 2h p Vr (4) 

where Vi = |iir| is the radial velocity of the mass-flow, 
E ~ pH is the surface density, and p is the mid-plane density 
of the disk. Hydrostatic equilibrium in the direction perpen¬ 
dicular to the disk mid-plane gives h ~ Cs/Hk where Cs is 
the local sound speed and VK(r) = r IlK(r') is the Kepler 
speed of the disk at radius r. Using the a prescription of the 
kinematic viscosity, v = a Cs h (Shakura & Sunyaev 1973), 
we have 

V h , , 

Ur ~ = a Cs — = a ^ Vk (5) 

which shows that vk ~ 10® Vr with h/r ~ 10~^, and a ~ 
0.1. This shows how slowly the matter flows in the radial 
direction relative to the Keplerian orbital speed. 

In the boundary region, the timescale of interaction be¬ 
tween the inner disk matter and the field lines tint — (D* — 
DK(nn))“^ = (cj - 1)“^ flK(nn)“^ where CU = fI*/DK(nn) 
is the fastness parameter. The escape speed at a given ra¬ 
dius Vesc = V^vk. To expel the matter from the system, the 
speed of the field lines should be greater than Vesc. This is 
satisfied outside the radius where cu = II*/f2K(rin) = y/2. 
At r = Cco the field lines and the matter move with the 
same angular speed. Since Dk oc , the propeller mech¬ 

anism could be sustained with an inner disk radius greater 
than (■\/2)^/®rco — 1.26 Vco. In the boundary region, the 
magnetic torque acting on the matter should accelerate the 
matter to the speed of the field lines within the interac¬ 
tion timescale, Lnt.This means that the magnetic stresses 
should accelerate the matter with the rate 17 ~ A17/tint 
where A17 — Q, — 17k (cin), which requires an angular mo¬ 
mentum transfer to the unit volume of the matter with the 
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rate 

2 A 2 
pr it pr 


An 


(n* — JIk) ^ 


= (a; - 1)^ pr^ r2K(nn)^. (6) 


Then, the magnetic torque per unit volume of matter, tb, 
should satisfy the critical condition 

T-B = ^ ~ (cu - if pr^ nKfinf (7) 

for a steady propeller effect. Apart from the r/h term on 
the left hand side, equation (7) looks similar to equation (2) 
which gives the conventional Alfven radius with the free-fall 
speed Vs ^ rflK- Nevertheless, for a given accretion rate at 
a radius r, the mass density in the disk is {r/h){vK/vi) ~ 
a~^{r/hf ~ 10 ^ times greater than the mass density in the 
spherical accretion. This means that the value at the 
conventional Alfven radius is about 5 orders of magnitude 
lower than that required for the propeller effect. To put it 
in other words, the inner disk radius in a steady propeller 
phase is likely to be much smaller than the Alfven radius 
given by equation ( 2 ). 

Using equations (3 - 5) and (7), we obtain the maximum 
inner disk radius in terms of the conventional Alfven radius 

2/7 


(cu- !)■ 




ta- 


( 8 ) 


Using the notation a_i = a/0.1, and (/i/r )_2 = (/i/r)/10 
equation ( 8 ) can be rewritten as 


Tin 3.8 X 10 


-2 2/7 

/A/ ' 


-V -1) 


-4/7 


4/7 




ta 


( 9 ) 


where hjr should be calculated at r = rin for the unper¬ 
turbed disk. Note that hjr oc has a weak de¬ 

pendence on both r and Min (see e.g. Frank et al. 2002). 
We would like to stress that equation (9) is valid only in 
the propeller phase, that is, when cj > %/2. In this regime, 
the rin/rA ratio increases with decreasing lo, and reaches a 
critical value corresponding to propeller-accretion transition 
when uj ~ y/2. , Using equation (3), we obtain 


~ 6.3 X 10 


-2 2/7 




( 10 ) 


For large ranges of Mn (10^^ — 10^® g s“^) and r (10^ — 10^° 
cm), the h/r ratio remains between ~ 2 x 10“® — 10“^. For 
Q = 0.1, Ta is found to be ~ 15 — 40 times greater than rin 
during the accretion-propeller transition. 

The boundary layer is a region of continual opening and 
reconnections of the field lines. From this region the matter 
leaves the system along the open field lines (see e.g. Lovelace 
at al. 1999 for a detailed description of the model). The outer 
boundary has a surface area A = 27rrin2/i. The mass per 
unit volume passing through this surface will rapidly gain 
an angular momentum ~ rBfint through interaction with the 
field lines within tint- The disk matter enters the boundary 
with radial speed Vi much smaller than its azimuthal speed 
t >0 ~ vk- Inside the boundary, that matter is accelerated 
in the azimuthal direction to the speed rinfl* = turinflK 
within tint = {ui — 1)”^ flK(nn)~^. Using equation (4) with 
Min = Mont, the total angular momentum transfer rate from 
the field to the outflowing matter becomes 


L — Tb tint Vr A — Afin (aj 1) Tin flK(Tin) 


( 11 ) 


when Tm — Tin > 1.26 Tco- This could be taken as a lower 
limit for the torque acting on the star, since we have ne¬ 
glected the effects of possible outflows returning back to disk 
which might lead to a pile up at the inner boundary. All 
these effects tend to further decrease the inner disk radius 
in the propeller phase. The inner disk radius in a steady- 
state cannot be built up at radii greater than Tin given by 
equation (9). This is because the required angular momen¬ 
tum for the matter to escape should be transferred within 
the interaction timescale. In particular, at the Alfven ra¬ 
dius, this needs a time interval that is about 10 ® interaction 
timescale. Initially, the matter could be stopped at ~ ta, 
but since the matter cannot be accelerated to high speeds 
there, the matter piles up at the inner disk. This extends the 
inner disk inwards until the critical condition for a steady 
propeller effect is satisfied at the inner disk radius. When the 
system reaches the steady-state, the disk regions outside the 
boundary region remain decoupled from the field lines, and 
thereby in an unperturbed steady thin disk state. We have 
proposed that this steady state could be achieved with an 
inner disk radius less than around Tin given in equation (9). 


2.3 The critical accretion rate for 
accretion-propeller transition 

The critical mass-flow rate of the disk, Min, for transition 
between the accretion and the propeller phases could be es¬ 
timated by setting Tin = 1.26 tco and uj = y/2 in equation 
(9). Using also equation (3), and Tco = we 

find 

Merit 4 X 10^" gs-^ a_i PjJ/® ^26 (12) 

where p 26 is the magnetic dipole moment of the star in units 
of 10^® G cm®, and P _3 is the spin period of the star in 
units of milliseconds. The propeller phase could prevail up 
to a maximum mass inflow rate of the disk given by equation 
( 12 ), which corresponds to critical accretion and disk lumi¬ 
nosities Lx = GMM/R and La = GMM/rin respectively, 
with M = Merit- 

Recently, Joadand et al. (2016) showed through detailed 
timing analysis that J1023-1-0038 slows down in the X-ray 
pulsar phase with a rate z> = —3.0413 x 10“®® Hz s“®) which 
is about 27% higher than the rate in the radio-pulsar phase 
{v = —2.3985 X 10“®® Hz s“®). This small fraction of increase 
in the torque when the system undergoes a transition from 
the radio pulsar to X-ray pulsar state might indicate that the 
magnetic dipole torque dominates the propeller torque in the 
radio pulsar phase. For a purely magnetic dipole torque, the 
dipole field strength on the pole of the star can be estimated 
through P - 6.4x 10®® G. For J1023-b0038, which has 
very similar properties to other transient sources, observed 
period and period derivative give B — 2 x 10® G (half of 
this at the equator) which is a typical field also for the other 
millisecond pulsars with shortest periods (see e.g. Lyne and 
Graham-Smith 2007). 

For a MSP with P = 1.69 ms and P ~ 2 x 10® G on the 
surface of the star, the condition ta = 1.26 Tco, is obtained 
with an accretion rate of ~ 6 x 10®® g s“® (equation (3)), 
which is about 4 orders of magnitude greater than the criti¬ 
cal accretion rate indicated by equation (12). It is seen that 
Merit oc This implies that the systems with smaller 
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periods have higher critical accretion luminosities. That is, 
for a given field strength, the systems with relatively low pe¬ 
riods undergo accretion-propeller transition at higher X-ray 
luminosities. This could be the reason for that the observed 
transitional MSPs to be short-period systems. The critical 
accretion rate given in equation (12) is consistent with the 
properties of MSPs during the transitions between the X-ray 
pulsar and the radio pulsar states (see Section 3 for further 
discussion). 


3 DISCUSSION AND CONCLUSIONS 

We have investigated the required conditions for a steady- 
state propeller phase, for which all the matter arriving at 
the inner disk radius rin > 1.26 rco can be expelled from the 
system with Vout > Vesc and Min = Mont- We have found 
that the propeller state could be established with an inner 
disk radius more than ~ 15 times smaller than the conven¬ 
tional Alfven radius depending on the dipole field strength 
and the period of the star, and the mass-flow rate of the 
disk. 

We derive rin with the assumption that the disk out¬ 
side the boundary has the properties similar to a standard 
thin disk. For greater surface densities due to some weak 
interactions between the held lines and the disk outside the 
boundary region, the inner disk would then settle down with 
an inner radius even smaller than the radius given by equa¬ 
tion (9) . We take the maximum possible magnetic stresses in 
our calculations; similarly, weaker magnetic stresses would 
give relatively small inner disk radii. Nevertheless, consid¬ 
ering very sharp radius dependence of the magnetic torques 
the actual inner disk radius, in the propeller regime, is likely 
to be near the maximum possible inner disk radius rin given 
by equation (9). 

When the orbital kinetic energy difference of the matter 
with vk and rinfl* is larger than the thermal energy density, 
Kelvin-Helmholtz instabilities become important for the gas 
to penetrate into the boundary region (Spruit and Taam 
1993) . The diffusion timescale of the disk into the boundary 
region tdiff ~ A7?/udifT where Udiff is the minimum of Av ~ 
rin(fl* — Hk) and Cs. For radii greater than 1.26rco, Av 
is much greater than Cs, and Vr <C Cs for a thin disk. This 
means that the inner disk can penetrate the boundary region 
without forming a pile-up in front of the boundary in the 
propeller phase. Furthermore, all the matter flowing into the 
boundary can be expelled from the boundary with the same 
rate as the rate of mass-flow into the boundary. This picture 
is consistent with our standard thin-disk assumption. 

In the numerical simulations, the efficiency and the 
time dependence of the outflows depend on some model 
assumptions (see Uzdensky 2004 for a review). In the 
present work, to estimate the critical condition for the pro¬ 
peller phase, we have adopted the basic principles of the 
disk-magnetosphere interaction model proposed initially by 
Lovelace et al. (1995), and developed further by Lovelace 
et al. (1999) to study the outflows in different regimes. The 
details of the propeller phase in the same model was studied 
through numerical simulations by Romanova et al. (2005) 
and Ustyugova et al. (2006). In these simulations, the mag¬ 
netosphere is dehned as the region in which all the matter 
and the closed held lines rotate together with the angular 


speed of the star while any plasma Hows are along the Held 
lines. The magnetospheric radius is taken to be equal to the 
radius at which ~ P + pv^ in the boundary where 

the thermal pressure, P, is negligible under the conditions 
we consider here. The matter inside the boundary is ac¬ 
celerated from quasi-Keplerian speeds. The inner disk ra¬ 
dius Tin, which is not known a priori, is determined from 
the simulations. The inner disk radius is always found to 
be close to the magnetospheric radius (rin — Vm) within a 
factor of ~ 3 in the propeller phase even when the systems 
show diHerent outHow morphologies (Ustyugova et al. 2006). 
In line with these results, the inner disk radius, rin, given 
in equation (7) remains close to r-m found by Lovelace et 
al. (1999) for all reasonable values of uj, except for when 
the system is very close to the propeller-accretion transition 
(rin = 1.26rco). From equations (7) and (8), it can be shown 
that Tin ~ (a; — l)“'‘'^^(r//i)^'^^rni. For r/h ~ 10^, the rin/rm 
ratio changes from 3.7 to 1 for the values of (cu — 1) from 1 
to 10. Nevertheless, for uj values greater than ~ 10, our rin 
decreases slightly below rm found by Lovelace et al. (1999). 
For instance, ui = 20 gives rin cs 0.7rm. 

The best sources to test our results are the recently dis¬ 
covered and clearly observed transitional MSPs which show 
transitions between the accretion powered X-ray pulsar and 
the rotational powered radio pulsar states. Two of these 
sources, namely PSR J1023-I-0038 (Archibald et al. 2015) 
and XSS J12270—4859 (Papitto et al. 2015), are observed to 
show X-ray pulsations in the sub-luminous X-ray states with 
Lx « a few 10®^ — 10®^ erg s“^. Properties of these X-ray 
pulsations strongly indicate that they are powered by mass- 
How onto the neutron star. These sources are observed in 
the radio pulsar phases when Lx ~ a few 10®^ erg s“^ with 
no X-ray pulsations. It seems that the radio pulsar phase is 
switched on at a disk mass-flow rate lower than ~ 10^® g 
s~^, and above this X-ray luminosity the sources are in the 
accretion phase. Since the accretion switches off the pulsed 
radio emission, the sources are very likely to be in the strong 
propeller phase when they emit radio pulses. Using the crit¬ 
ical condition rin — 1.26 Too for accretion-propeller transi¬ 
tion, we can estimate the X-ray luminosity originating from 
the accretion onto the star. Lx, just before the transition 
to the propeller state and the X-ray luminosity of the disk, 
Ld, just after the transition. The period P = 1.69 ms of the 
transitional MSPs corresponds to rco — 2.4 x 10® cm. For a 
neutron star with radius R = 10® cm and mass M = 1.4M0, 
the transition luminosity Lx = GMM/R ~ 1 x 10®® erg s“® 
corresponds to an accretion rate M ~ 5 x 10®^ g Imme¬ 
diately after the transition to the propeller phase, the X-ray 
luminosity decreases to Ld = GMMin/2rin ~ 1.5 x 10®^ erg 
s“®. Note that these estimates are independent of the model 
details, except for the assumed inner disk position during 
the transition. 

Is Tin indeed close to 1.26 rco for M ~ 5 x 10®^ g s“®? In 
other words, is the critical accretion rate corresponding to 
this rin in the models, consistent with this critical accretion 
rate? In our model, we can check this from equation (12). 
For P = 1.69 ms and B ^ 10® — 10® G, the values of the 
critical accretion rate Mcut calculated from equation (12) 
are found to be consistent with the transition rate indicated 
by the observations. Considering the possibilities of the in¬ 
teraction between the pulsar wind and the inner disk, and 
the magnetic heating in the disk-fleld interaction region, the 
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observed X-ray luminosity in the radio pulsar phase could 
be taken as an upper limit to the standard disk luminosity. 

Our result for the inner disk radius is applicable in the 
propeller phase (r^o < 1.26rin). In the spin-down with accre¬ 
tion phase, when r^o < rin < 1.26 Cco, the matter cannot be 
thrown out with sufficiently high speeds and piles up at the 
inner disk. This could lead to repetitive switches between the 
accretion and the propeller phases on dynamical and viscous 
timescales of the innermost disk. These transitions might be 
responsible for the observed rapid recurrent changes between 
two different X-ray modes of the transitional MSPs with lu¬ 
minosity ratios of ~ 5 — 7 (see Linares 2014 for a review). It 
is remarkable that Lx/Ld — 6 when the inner disk radius is 
Tin ~ 1.26 Cco — 3 X 10® cm. This phase, neglecting the pro¬ 
peller effect, was studied in detail by D’Angelo and Spruit 
(2012). In this regime, the total torque acting on the star 
is also expected to be rather different from that in the pro¬ 
peller phase. Because, the pile up at the inner disk required 
to satisfy rin = r^o increases the viscous stresses consid¬ 
erably. Furthermore, the interaction timescale of the inner 
disk matter and the field lines increases as rin approaches to 
rco, which makes the transition to the spin-up phase rather 
complicated together with uncertainties in the radial width 
of the boundary region. We will investigate the properties 
of the spin-down phase with accretion and the transition to 
the spin-up phase in an independent work. 

In summary, we have estimated the maximum inner disk 
radius in the propeller phase and the critical accretion rate 
for the accretion-propeller transition for magnetized stars 
with optically thick geometrically thin standard accretion 
disks. We have shown that a steady propeller phase could 
be achieved with an inner disk radius much lower than the 
conventional Alfven radius. In this phase, the inner disk ra¬ 
dius decreases with increasing mass-flow rate of the disk. 
At a critical accretion rate, corresponding to rin — 1.26 rco, 
the system makes a transition from the propeller to the ac¬ 
cretion phase. This accretion rate is much lower than the 
rate that can balance the magnetic and the viscous stresses 
close to the co-rotation radius. The transition from the ac¬ 
cretion into the propeller phase is likely to be accompanied 
by the switch-on of the pulsed radio emission, while all the 
mass flowing into the boundary region is being effectively 
thrown out along the open field lines. Our results, that are 
in agreement with the properties of transitional MSPs, can 
also be extended to other magnetized stars accreting from 
geometrically thin accretion disks. 
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